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Abstract. We consider an open quantum Fermi-system which consists of a single 
degenerate level with pairing interactions embedded into a superconducting bath. The 
time evolution of the reduced density matrix for the system is given by Linblad master 
equation, where the dissipators describe exchange of Bogoliubov quasiparticles with 
the bath. We obtain fixed points of the time evolution equation for the covariance 
matrix and study their stability by analyzing full dynamics of the order parameter. 
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The Bardeen, Cooper, and Schrieffer (BCS) theory of superconductivity pQ is 
based on the simple Hamiltonian but it captures the essential physics not only for 
superconductivity of electrons in metals, but also the superconductivity of atomic nuclei, 
nuclear matter, neutron stars [2], and cold atomic Fermi-gases [3]. Here, based on the 
Lindblad master equation, we propose the extension of the BCS theory to the open 
quantum system. We determine the order parameter - the density of Cooper pairs, 
and optionally the order parameters of a superconducting reservoir, self consistently. 
A simple model of a single spinful fermionic level embedded into a fermionic reservoir 
is proposed. For a self consistent treatment of the reservoir we recover the standard 
mean-field superconducting phase transition in the grand-canonical (equilibrium) state, 
whereas for a fixed state of the reservoir, we find that the state of the system follows 
the state of the reservoir, being either superconducting or normal. 

We consider a Lindblad master equation for the time-dependent reduced density 
matrix p(t) [I] 

d M 

£ = -i[H, p} + J2 (2£,P4 " {LlL v , p}) , (1) 

u=l 

where Hamiltonian H is given as 

H = e(aja-)- + a\a±) + Uala^a^ . (2) 

The Hamiltonian consists of a single degenerate level with BCS type pairing interaction; 
al,a a are standard fermionic creation/annihilation operators (er =f, 4), and L M are 
some Lindblad operators which will be specified later. Introducing (time-dependent) 
averages (A) = txAp(t), and using Wick theorem we approximate ^ by the mean-field 
Bogoliubov-de Gennes Hamiltonian 

H = e(ata-r + ajajj + A(e lx a±a^ + e~ lx a|a|), (3) 

where 

Ae ix = U(a\a\) (4) 

is the complex order parameter. 

The Hamiltonian ^ can be diagonalized by the canonical Bogoliubov 
transformation 

«i = — e lx cos(0)a-|- + sin(0)a|, (5) 



a 



2 = e 1X cos(0)a| + sin(0)a-f-, (6) 

where tan (20) = — — , and cti are Bogoliubov quasiparticles, which satisfy standard 
anticommutation relations. The diagonal form of the Hamiltonian is 

H = Ve 2 + A 2 + a\a 2 ). (7) 

We assume that our system is embedded into a superconducting bath described by a 
macroscopic gas of Bogoliubov quasiparticles, which can be exchanged with the system. 
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This process is quite generally described by a combination of the following M 
Lindblad operators 

-e ir > cos(6 l )a t + sin(6»)a| 
-e~ iri cos(8)a\ + sin(6>)a^ 



Li 

u 



r 2 



(8) 



Ti (e ir? cos(fl)a 4 + sin(#)a} 
L , = \/V 2 { <r ir? cos(#)a{ + sin(0)a t 



Superconducting properties of the bath are not necessarily the same as of the embedded 
system, so the angles 9, and r], may be different from 0, and Xi respectively. The coupling 
constants Tip are connected to the Fermi-Dirac distribution of the normal modes of the 
Bogoliubov-de Gennes Hamiltonian ([3]): 

rWW), r 2 = 7/ , f = i + J^ , (9) 

where 7 is a parameter controlling the strength of the system-bath coupling. 

It is convenient to rewrite the Lindblad equation in terms of Hermitian Majorana 
fermions 



w\ = + aj, W2 
satisfying {wj,w k } = 25 j>k , j,k = 



a\), 



(10) 



1, ... ,4. Since our master equation is quadratic in 
terms of Wj, we obtain a closed set of equations for the covariance matrix (wjWk) = 
8j,k — iZj jk (t): 

dZ 

~dt 

Here Z(t) is a real, anti-symmetric 4x4 matrix, X and Y are real matrices 



-X T Z - ZX + Y. 



(11) 



X 



Y = 2 7 (1 - 2f) 



( 7 

e 

A sin% 
V — A cosx 

/ 

-cos 29 

sin r) sin 29 
\ — cos 7] sin 29 



7 



— e — Asinx Acosx\ 
7 A cos x A sin x 

-Acosx 7 — e 

-A sin x e 
cos 29 


— cos 1] sin 29 
sin 7] sin 29 — cos 29 



(12) 



— sin i] sin 29 cos r\ sin 29 \ 
cos 77 sin 2$ sin 77 sin 2$ 



cos 20 




(13) 



/ 



Note that eq. (11) is non-linear, as A and x depend again on covariances through the 
relation Q. If in addition, we determine the Lindblad operators self-consistently, which 
means that the bath has the same properties as the embedded system, then we should 
also set tan 29 — ——, (9 = 0), and rj = x- 

Fixed points of the flow (11), which are solutions of the continuous Lyapunov 
equation 



X T Z + ZX = Y, 



(14) 
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Figure 1. The phase diagram of the magnitude of the order parameter A versus the 
inverse temperature. A > signals the superconducting phase. Blue curve corresponds 
to self-consistently determined bath, whereas black and red curves correspond to fixed 
superconducting baths, with Abath = — e tan 26 (indicated by horizontal dashed lines), 
with 9 = — 7r/20 and 8 — —tt/7, respectively. Other parameters are e = 1, 7 = 1, U — 
-3. 



determine the stationary states of the system. These stationary states may not be 
unique because of the non-linearity. However, one can show that a stable fixed point is 
unique. It follows from the fact that all linear relaxation rates, i.e. eigenvalues of the 
matrix X (12), £1,2,3,4 = 7 ± i^A 2 + e 2 , have strictly positive real parts for 7 > 0. One 
can further show that all solutions of Eq. ( 14 ) are of the form 



(15) 



which is specified by only three real variables zi, z%, Z3. The Lyapunov equation should 
be solved self-consistently (J4J), i.e. 

2A . 2A 



/ 





Zl 


Z 2 


zs \ 









Z3 


-zi 




-z 2 


-z 3 





Zl 


V 


-z 3 


Z-2 


-zi 


/ 



Z2 = -— smx, 



z 3 



U 



cosx- 



(16) 



Let us now consider two possible cases of the bath, i.e. two possible choices of angles 

(i) Fixed bath. If the bath is considered fixed, we can set r\ := without loss of generality, 
so the equation (16) together with the Lyapunov equation (14) results in the conditions 

(7 sin x 



e cosx) sin 26> = A cos 26* , (17) 
2f) ( ( ( 7 2 + A 2 ) cos x + 7e sin X ) sin 29 - Ae cos 26) = 2 A ( 7 2 + A 2 + e 2 ) . 



For normal, non-superconducting bath, sin 29 = 0, and there exists only a trivial solution 
A = for the system. For a superconducting bath, sin 29 7^ 0, the system is always 
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superconducting as well, for all inverse temperatures /3, as the only solution (18) has 
A^O. In other words, for a superconducting bath the embedded system cannot be in 
the normal state no matter how small the parameter 7 is. 

(ii) Self- consistent bath. The bath can be considered self-consistent in the following way. 
We associate the Lindblad operators with the normal mode creation and annihilation 
operators such that 9 = 0, 77 = %■ Inserting these assumptions in the Lyapunov equation 
(14) we get 



z\ 



z 2 



Z3 



etanhQ/VA 2 + e 2 ) 
VA 2 + e 2 



A sin(x) tanh Q/VA 2 + e 2 ) 
VA 2 + e 2 



A cos(x) tanh Q/VA 2 + e 2 ) 
VA 2 + e 2 



The consistency conditions (16) now result in 



£/tanhQ/VA 2 + e 2 ) 
2VA 2 + e 2 

We always have the trivial solution A 
where 



+ 1 A = 



(19) 



P c = — artanh 
e 



0, however for temperatures smaller than 1//3 C 



(20) 



we also find superconducting solution with A / 0. Obviously, critical temperature 
only exists for —U > 2e. Eqs. (19, 20) are exactly the result that we get in standard 
grand canonical ensemble. The bath-self-consistency conditions completely eliminate 
dependence of the properties of the stationary state on the system-bath coupling 7. 

In figure [l] we plot a phase diagram A(/3), for both cases (i,ii). Note that only if 
Abath := —e tan 20 is smaller than J jU 2 — e 2 then we may have A = Abath for some 
temperature 1/(3. 

Let us now 
We shall 



So far we have investigated fixed points of the non-linear flow (11). 



address the question of their stability by investigating the full dynamics, 
only focus on the case (ii) of self-consitently determined baths, which possesses non- 
unique stationary solution below the critical temperature. The dynamical equations for 
Zj(t),j = 1,2,3 follow directly from equation (11) with the ansatz (15). By means of 
self-consistency equations (16) we replace zzif), 23 (t) by A(t), x{P), resulting in a system 
of three non-linear differential equations 



dA 

d^ 

d^i 
dt" 

dx 
dt 



/f/tanh(|Ve 2 + A 2 ) 
27A %== ' 



2 7 



\ 2 V / ?TA 
tanh (|Ve 2 + A 2 ) 



+ 1 



v 7 ? 



+ A 2 



Uz x + 2e. 



(21) 

(22) 
(23) 
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Figure 2. Time-dependent order parameter A(£) as a function of time, for j3 < j3 c (red 
curves) and /3 > C (black curves), and for an attractive potential U < (full curves) 
and a repulsive potential U > (dashed curves). We plot some typical trajectories 
starting from green points and ending in one of the two stable fixed points (blue). 
Numerical values of parameters are e=l,7 = l,|C/|=3. 



Note that the first equation (21) is independent from the other two, and yields a closed 
first order differential equation for the magnitude of the order parameter A. Writing it 
as dA/dt = G(A) we can easily study the stability of two possible fixed points Ai, and 
A 2 (G(Aj) = 0, j = 1, 2, Ai = 0, and A 2 ^ if /3 > (3 C ). For the trivial fixed point we 
find 

G'(0) = -2 7 (— tanh ( ^ ) + 1^ , (24) 



2e V 2 , 

which means that the non-superconducting state is stable G'(0) < if (3 < f3 c and 
unstable G'(0) > if (3 > (3 C . As for the second, non-trivial fixed point A 2 , one can 
easily show that it is always stable G'(A 2 ) < if it exists, i.e. if /3 > (3 C . The flow of the 
order parameter A(£) in different cases is illustrated in Fig. [2] We note that the other 
two phase-space variables Zi(t),x(t) are completely enslaved by the order parameter 
A(t), so they again converge either to their trivial (non-superconducting) or non-trivial 
(superconducting) fixed point values. 

In conclusion, based on the Lindblad master equation we proposed the extension of 
the BCS theory to open Fermi systems exchanging Bogouliubov quasiparticles with the 
bath. We derived the equations of motion for the covariance matrix and found the fixed 
points of the flow equations. If the bath is considered self-consistently with the system, 
the results becomes equivalent to the grand canonical ensemble and all dependences on 
the Lindblad dissipators are eliminated. If the superconductivity of the bath is fixed, 
the system remains in the superconducting state for all values of temperature no matter 
how small the system-bath coupling is. We performed the stability analysis of the fixed 
points and found that below the critical temperature the fixed point which corresponds 
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to the normal phase is not stable, whereas the superconducting solution is stable. Note 
that our results are closely related to an exact treatment of open BCS model in quasi- 
spin formulation [5]. However, in Ref. [5] the (quasi-)particles cannot be exchanged with 
the system, so the study of thermodynamic limit is more subtle. We thus believe that 
our single-level formulation provides a minimal model of open BCS quantum dynamics 
and should serve as the first step in approaching the non-equilibrium open BCS models 
with several different temperature/particle reservoirs. 
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